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We accept the implicit challenge of A. Uhlmann in his 1994 paper, “Parallel Lifts and Holonomy 
along Density Operators: Computable Examples Using 0(3)-0rbits,” by, in fact, computing the 
holonomy invariants for rotations of certain n- level Gibbsian density matrices (n = 2,..., 11). From 
these we derive, by the tracing operation, the associated geometric phases and visibilities, which 
we analyze and display. We then proceed analogously, implementing the alternative methodology 
presented by E. Sjoqvist et al in their letter, “Geometric Phases for Mixed States in Interferometry” 

(Phys. Rev. Lett. 85, 2845 [2000]). For the Uhlmann case, we are able to also compute several 
higher-order holonomy invariants. We compare the various geometric phases and visibilities for 
different values of n, and also directly compare the two forms of analysis. By setting one parameter 
(a) in the Uhlmann analysis to zero, we find that the so-reduced form of the first-order holonomy 
invariant is simply equal to ( —l) n+1 times the holonomy invariant in the Sjoqvist et al method. 

Additional phenomena of interest are reported too. 
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X Summary 


I. INTRODUCTION 

In a recent paper, Sjoqvist et al [Q provided “a new formalism of the geometric phase for mixed states in the 
experimental context of quantum interferometry”. Only in passing, did these authors note that “Uhlmann was 
probably the first to address the issue of mixed state holonomy, but as a purely mathematical problem”. 

Interested in possible (previously uninvestigated) relationships between the work of Sjoqvist et a/[lj and that of 
Uhlmann |^||, the present author compared the two approaches in terms of two spin-^ scenarios |i ] . In the first 
of these, the spin-1 systems undergo unitary evolution along geodesic triangles Q], while in the second the unitary 
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evolution takes place along circular paths ||. In ||, Uhlmann had also proposed an additional scenario, which was 
(initially) left unanalyzed in M. It involves “the Gibbsian states of the form”, 


„a n 


7 


p = 


trace e “ n 


T 


(i) 


which fill for a given value of a a 2-spliere called S? if ft runs through all directions in 3-space. Now, for the starting 
point of the unitary evolution, one sets it = (0, 0,1), while rf = (0, sind, cos 9) is chosen as rotational axis. The curve 
of state evolution is given by 


</> -> U(4>)poU{-4>), 17(0) = e-^( sinej *+ cosej *), 

and the associated parallel lift of this curve with initial value p 1 / 2 is 

0 - f/(0)pJ / V(0), V(0) = e** 6 , 

where 

1 


H = cos 9J Z + a sin 9 J y , and a = 


cosh ■ 


( 2 ) 


(3) 


(4) 


It is possible to regard V (0) as a rotation with angle 


and rotation axis 


= K(f>, k = 7 cos 2 9 + a 2 sin 2 9 < 1 


—> . sind a cos 9. 

S • )• 

hi K 


The holonomy invariant can then be written as 

(-l 2j )Po /2 e 2 "^Po /2 = (-1 2j )pl /2 e 2niK ^ p 1 ^ 2 . 


(5) 


( 6 ) 


(7) 


(All the preceding equations are directly adopted from ||.) 

In this paper, to begin, we compute the trace of this invariant (7!) for all j = J, 1, |,..., |. We plot the arguments 
of these traces, that is the corresponding geometric phases (lj), in Figs, [j] - |l0[ (sec. ||), and their absolute values, 
that is the visibilities ( Vj ), in Figs. (sec. p!| ). (For j = |, the results are equivalent to the first (non-Gibbsian) 

set of models in which was compared with the analyses of Sjoqvist et al Q in [Q.) Let us note that the value 
a = 0 corresponds to the fully mixed (classical) state, while a = ±oo correspond to pure states. 


In sec. IV, we compare certain aspects of these results across the ten distinct values of n. All the various results up 
to this point can be co nsider ed as first-order in nature. In sec. |v|, on the other hand, we compute certain higher-order 
invariants. In secs. VI - VIII , we perform analyses precisely analogous to those in secs. p| - ^v| , but now in terms not 
of the approach of Uhlmann to mixed state holonomy, but that of Sjoqvist e£ aZ [jl| . In our final analytical section 
(sec. IX), before the summary (sec. [x]), we directly compare results obtained by the two different procedures. 

One of our main findings is that if one sets the implicit parameter a in the Uhlmann holonomy invariant ( 0 ) to 
zero, then its resultant trace is simply equal to (—l) n+1 times the holonomy invariant [0, eq. (15)] yielded by the 
methodology of Sjoqvist et al. 
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II. UHLMANN GEOMETRIC PHASES FOR GIBBSIAN TV-LEVEL SYSTEMS (TV = 2,, 11) 



FIG. 1. Uhlmann geometric phase for Gibbsian spin-^ systems 
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FIG. 2. Uhlmann geometric phase for Gibbsian spin-1 systems 
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FIG. 3. Uhlmann geometric phase for Gibbsian spin-^ systems 
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FIG. 4. Uhlmann geometric phase for Gibbsian spin-2 systems 



FIG. 5. Uhlmann geometric phase for Gibbsian spin-§ systems 
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FIG. 6. Uhlmann geometric phase for Gibbsian spin-3 systems 



FIG. 7. Uhlmann geometric phase for Gibbsian spin-4 systems 
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FIG. 8. Uhlmann geometric phase for Gibbsian spin-4 systems 



FIG. 9. Uhlmann geometric phase for Gibbsian spin-| systems 
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FIG. 10. Uhlmann geometric phase for Gibbsian spin-5 systems 


III. UHLMANN VISIBILITIES FOR GIBBSIAN A-LEVEL SYSTEMS (A = 2,..., 11) 



FIG. 11. Uhlmann visibility for Gibbsian spin-i systems 



















FIG. 12. Uhlmann visibility for Gibbsian spin-1 systems 



FIG. 13. Uhlmann visibility for Gibbsian spin-| systems 
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FIG. 14. Uhlmann visibility for Gibbsian spin-2 systems 



FIG. 15. Uhlmann visibility for Gibbsian spin-| systems 
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FIG. 16. Uhlmann visibility for Gibbsian spin-3 systems 



FIG. 17. Uhlmann visibility for Gibbsian spin-| systems 
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FIG. 18. Uhlmann visibility for Gibbsian spin-4 systems 
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FIG. 20. Uhlmann visibility for Gibbsian spin-5 systems 

IV. COMPARISONS ACROSS GIBBSIAN A-LEVEL SYSTEMS OF UHLMANN GEOMETRIC PHASES 

AND VISIBILITIES 


In the spin-1 case (Figs. |l] and 0 ) , we have that the geometric phase 71 is the arctangent of ^, taking into account 
which quadrant the point (x, y) lies in, where 


x = — coshi7TK, 

Also, the visibility in this spin-1 case is given by 


cos 9 sin 7 tk tanh % 

y = - 1 


Vl 

2 



4 sinh 2 inn 

c 


( 8 ) 


(9) 


where 


C = 3 — cos 29 + 2 cos 2 6 cosh a, (10) 

and k is as defined in <D- The comparable results in the spin-1 case (Figs. and [l2|) are similarly based on 


x = 


: (2{e a + 2e 21 ™ + e Q+4i7rK ) cos 2 9 + (2 + 2e 4i ™ + e“(l + e 2l7r,t ) 2 )sech 2 f sin 2 

4k 2 (1 + 2 cosh a) 


and 


V = 


e“ cos 9 sin 27 tk 
k + 2kcosq: 


( 11 ) 


( 12 ) 


One feature distinguishing the (strikingly similar) visibility plots (Figs. |lT] - 20) from one another is that as j increases, 
Vj decreases at the (boundary) point a = 5, 9 = For instance this value is .871618 in Fig. |TT| and .498776 in Fig. [Ii]. 
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In Fig. |2l], we plot 7 j (j = ...,) versus 9 , holding a fixed at 1. Curves for lower-dimensional Gibbsian systems 
here strictly dominate those for higher-dimensional systems. (As a increases, however, above certain thresholds for 
6, this simple monotonic behavior vanishes [cf. Fig. p4j|.) 



FIG. 21. Uhlmann geometric phases for n-level Gibbsian systems (n = 2,..., 11) holding a = 1. The curve for n — 2 
dominates that for n = 3, which dominates that for n = 4,... 


In Fig. ^2], we 
the ten curves is 


“reverse” this scenario, now holding 9 fixed at and letting a vary over [0,5]. The monotonicity of 
completely analogous to that in Fig. with curves for lower j dominating those for higher j. 



FIG. 22. Uhlmann geometric phases for n-level Gibbsian systems (n = 2,... , 11) holding 9 fixed at The curve for n = 2 
dominates that for n = 3, which dominates that for n = 4,... 

In Fig. we hold cr at 2 and plot the visibilities for the ten spin scenarios. Again, as in Figs. |2l| and ^2[ curves 
for lower values of j dominate those for higher values. 
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FIG. 23. Uhlmann visibilities for n-level Gibbsian systems (n = 2,..., 11) holding a fixed at 2. The curve for n = 2 
dominates that for n = 3, which dominates that for n = 4,... 

While in Fig. |2l] the inverse temperature parameter a was fixed at 1, in Fig. ^4], it is held at 2. Now the same 
simple monotonicity with j observed in the preceding three figures, holds below 9 « .65. However, it is lost at higher 
values of 9 , with the curves for the higher spin states oscillating from -f to ^ at lower values of 9 than do the curves 
for some of the lower spin states. 


V 



FIG. 24. Uhlmann geometric phases for n-level Gibbsian systems (n = 2,..., 11) holding a = 2. Only below 6 « .65 does 
the curve for n = 2 dominate that for n = 3, which dominates that for n = 4,... 

In Fig. H||, in which 9 is fixed at similar behavior to that observed in Fig. ^ takes place. Below a ~ 1.43, the 
simple monotonicity with j holds, then the curves begin to jump from - j to f with the curves for higher j jumping 
first. 
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FIG. 25. Uhlmann geometric phases for n-level Gibbsian systems (n = 2,..., 11) holding 8 = j. Only below a « 2 does the 
curve for n = 2 dominate that for n = 3, which dominates that for n = 4,... 

In Fig. ^6], holding 9 = we plot the various visibilities over the range a € [1,4]. At a = 2.2, the values of Vj 
monotonically decline from j = \ to |. 



monotonically decline as n increases. The highest peak belongs to n = 2. 


V. HIGHER-ORDER UHLMANN HOLONOMY INVARIANTS 


A. Traces of powers of Uhlmann holonomy invariants 


All the analyses above have been conducted on the basis of the trace of the holonomy invariant ([^). In addition, 
however, the traces of the higher powers of ( 0 ) are also invariants || eq. ( 102 )] (all reducing, however, for pure states, 
to simple powers of the “Berry phase” ||, eq- (112.5)]) (cf. [||). In Fig. [27], we show the argument of the second power 
of the holonomy invariant for the spin-i case (n = 2 ), which we express in the form of the arctangent of -, taking 
into account which quadrant the point (x, y) lies in, where 


x = (C ( — 1 + (1 + cosh a) cosh 2 i 7 rK)sech 2 ^ — 8 sinh 2 intij , y 


2C 

The corresponding absolute value (Fig. 

, 9 „ 2cosh2i7TK 

x 2 + y 2 = l- — - 7 — 

1 + cosh a 


2in cos 9 sinh a sinh 2inn 

C ' 


is the square root of 


1 4 a 4(2 + cosh adsech 2 f sinh 2 inn 16 sinh 4 Am 

-sech 4 — +-^-+ C2 ■ 


(13) 


(14) 
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(In some of the later [computationally-intensive] figures in this section, it proved essentially necessary to omit the 
[otherwise obvious] axes labels, due to some quite distinct peculiarities of the graphics facilities for the local installation 
of MATHEMATICA on the more powerful workstations.) 




In Figs. and |3(i| are shown the analogous quantities for Gibbsian spin-1 systems. 
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FIG. 31. Argument of the trace of the third power of the holonomy invariant (|7|) for the f/iree-level Gibbsian systems (n = 3) 
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FIG. 33. Argument of the trace of the second power of the holonomy invariant (|ij) for the four -level Gibbsian systems (n = 4) 



FIG. 34. Absolute value of the trace of the second power of the holonomy invariant (Q) for the four -level Gibbsian systems 
(n = 4) 
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FIG. 35. Argument of the trace of the third power of the holonomy invariant (J7]) for the four-level Gibbsian systems (n = 4) 



FIG. 36. 
(n = 4) 



1.5 


3 0 


Absolute value of the trace of the third power of the holonomy invariant If for the four -level Gibbsian systems 
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FIG. 39. Argument of the trace of the second power of the holonomy invariant (]7j) for the five -level Gibbsian systems (n = 5) 
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FIG. 41. Argument of the trace of the third power of the holonomy invariant (^) for the five -level Gibbsian systems (n = 5) 


FIG. 42. 
(n = 5) 
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FIG. 43. Argument of the trace of the fourth power of the holonomy invariant (Q) for the five -level Gibbsian systems (n = 5) 
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FIG. 45. Argument of the trace of the fifth power of the holonomy invariant (H) for the five -level Gibbsian systems (n = 5) 
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FIG. 47. Argument of the trace of the second power of the holonomy invariant (|[7) for the sis-level Gibbsian systems (n — 6) 
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FIG. 49. Argument of the trace of the third power of the holonomy invariant (bj) for the sis-level Gibbsian systems (n = 6) 



FIG. 50. Absolute value of the trace of the third power of the holonomy invariant ([?]) for the sra-level Gibbsian systems 


























FIG. 51. Argument of the trace of the fourth power of the holonomy invariant (!7) for the six-level Gibbsian systems 


(n = 6) 
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FIG. 54. 
(n = 6) 



Absolute value of the trace of the fifth power of the holonomy invariant ([?]) for the 


SKC-level Gibbsian systems 
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B. Cross-comparisons of Uhlmann holonomy invariants 



FIG. 55. Ratio of the argument of the trace of the second power of the holonomy invariant 
of the first power for the eight-level Gibbsian systems 


to the argument of the trace 



FIG. 56. Ratio of the absolute value of the trace of the second power of the holonomy invariant 
the trace of the first power for the eight-level Gibbsian systems 


© 


to the absolute value of 
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FIG. 57. Ratio of the argument of the trace of the third power of the holonomy invariant (|?J) to the argument of the trace 
of the second power for the eight-level Gibbsian systems 



FIG. 58. Ratio of the absolute value of the trace of the third power of the holonomy invariant (7) to the absolute value of 
the trace of the second power for the eight-level Gibbsian systems 


C. Eigenvalues of Uhlmann holonomy invariant for n = 2,3,5 

To supplement our consideration above of the traces of various powers of the Uhlmann holonomy invariant (Jt|) , here 
we examine the eigenvalues themselves. (Both these traces and eigenvalues are themselves invariants, of course.) In 
Figs. |5^ and [io] are shown the arguments of the dominant and subordinate eigenvalues of the holonomy invariant (]7j) 
for the case n = 2, and in Figs. |bl] and ^2|, the corresponding absolute values. 
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FIG. 59. 
matrices 



Argument of the dominant eigenvalue of the Uhlmann holonomy invariant ([^) for the 


Arg[>] 


fuio-level Gibbsian density 



FIG. 60. Argument of the subordinate eigenvalue of the Uhlmann holonomy invariant (^) for the truo-level Gibbsian density 
matrices 
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In Figs. and we plot the arguments of the dominant and subordinate eigenvalues of the holonomy invariant 
( 0 ) for the case n = 3 and in Figs. and |0, the absolute values of the two corresponding eigenvalues. The argument 
of the intermediate (in absolute value) one of the three eigenvalues appears to always be identically zero, that is, 
this intermediate eigenvalue is real (Fig. |66| ). It also appears numerically, on the basis of further analyses we have 
conducted, that this result is generalizable to the proposition that for odd n, the central/most intermediate one (that 
is, the [^]-th) of the n eigenvalues is real. (Also, the centrally located diagonal entry of the holonomy invariant for 
odd n appears to be real, in general.) 
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FIG. 63. 
matrices 


FIG. 64. 
matrices 


Arg[A] 




Argument of the dominant eigenvalue of the Uhlmann holonomy invariant (f7|) for the fftree-level Gibbsian density 
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density matrices 



FIG. 66. Intermediate eigenvalue of the Uhlmann holonomy invariant ()7|) for the f/iree-level Gibbsian density matrices 
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FIG. 67. Absolute value of the subordinate eigenvalue of the Uhlmann holonomy invariant (|7|) for the three -level Gibbsian 
density matrices 

In Figs. |6^-[fl], we show the arguments of the eigenvalues of the Uhlmann holonomy invariant ( 0 ) in the case n = 5. 
(The case n = 4 proved more intractable in nature.) The third eigenvalue is simply real. In Figs. 0 - 00-0 are 
shown the corresponding absolute values. In Fig. 0 is shown the value of the (real) intermediate eigenvalue. 



FIG. 68. 
matrices 



Argument of the leading eigenvalue of the Uhlmann holonomy invariant rt7|) for the five -level Gibbsian density 
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FIG. 69. Argument of the second leading eigenvalue of the Uhlmann holonomy invariant 
matrices 


for the five-level Gibbsian density 



FIG. 70. Argument of the second smallest eigenvalue of the Uhlmann holonomy invariant (Q) for the five -level Gibbsian 
density matrices 


38 
















FIG. 71. Argument of the last eigenvalue of the Uhlmann holonomy invariant (pj) for the five-level Gibbsian density matrices 



FIG. 72. Absolute value of the leading eigenvalue of the Uhlmann holonomy invariant (|t|) for the five -level Gibbsian density 
matrices 
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FIG. 73. Absolute value of the second leading eigenvalue of the Uhlmann holonomy invariant 
density matrices 


1 


for the five-level Gibbsian 



FIG. 74. Central eigenvalue of the Uhlmann holonomy invariant (t|) for the five -level Gibbsian density matrices 
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for the five -level Gibbsian 


density matrices 



FIG. 76. Absolute value of the last eigenvalue of the Uhlmann holonomy invariant ( |t| I for the five -level Gibbsian density 
matrices 


VI. SJOQVIST ET AL GEOMETRIC PHASES FOR GIBBSIAN A-LEVEL SYSTEMS (A = 2,..., 11) 

In the approach of Sjoqvist et al |IJ every pure state that diagonalizes the initial density matrix is parallel transported 
separately. By weighting the associated holonomy of each such pure state by the corresponding eigenvalues, we obtain 
the measure proposed in Q “in the experimental context of quantum interferometry”. Since it is essentially trivial 
to obtain the pure states that diagonalize the initial (diagonal) n-level Gibbsian density matrix and the associated 
eigenvalues, and since Uhlmann has presented the holonomy invariant (the trace of which is the Berry phase) 

W 27r W 0 * = (-l) 2 Vj /2 e 2,ricoseJ Vo /2 = e - 2nirn{1 - cose) \m >< m\, (15) 
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for the circular evolution of pure states, we can directly compute the Sjoqvist et al geometric phase ( 7 ) and visibility 
(v) for the n-level Gibbsian density matrices. These are presented in Figs 0 to | 8 t] and Figs. |7| to [)(| The trace of 
the weighted (by the eigenvalues) sum of the holonomy invariants © for n = 2 is 


a 


cos (n cos 9) + i sin (n cos 9) tanh —, 


while for n = 3, it is 


1 + 2 cos (2-7T cos 9) cosh a + 2i sin (27 t cos 9) sinh a 
1 + 2 cosh a 


for n = 4, 


3i7r cos 6 


(i 


0 a.-\- c li'K cos 6 


)(i 


D 2o:+4z7r cos 6 \ 


(1 + e“)(l + e 2 “) 


(16) 


(17) 


(18) 


and for n = 5, 

1 + 2 cos (2-7T cos 9) cosh a + 2 cos (4 -tt cos 9) cosh 2 a + 2 i sin (27r cos 9) sinh a + 2i sin (4-7T cos 9) sinh 2a 

1 + 2 cosh a + 2 cosh 2a 


(19) 



FIG. 77. Sjoqvist et al geometric phase for Gibbsian spin-^ systems 
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FIG. 78. Sjoqvist et al geometric phase for Gibbsian spin-1 systems 
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FIG. 79. Sjoqvist et al geometric phase for Gibbsian spin-f systems 
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FIG. 80. Sjoqvist et al geometric phase for Gibbsian spin-2 systems 



FIG. 81. Sjoqvist et al geometric phase for Gibbsian spin-| systems 
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FIG. 82. Sjoqvist et al geometric phase for Gibbsian spin-3 systems 
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FIG. 83. Sjoqvist et al geometric phase for Gibbsian spin-| systems 
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FIG. 84. Sjoqvist et al geometric phase for Gibbsian spin-4 systems 
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FIG. 85. Sjoqvist et al geometric phase for Gibbsian spin-f systems 
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VII. SJOQVIST ET AL VISIBILITIES FOR GIBBSIAN TV-LEVEL SYSTEMS (V = 2,..., 11) 



FIG. 87. Sjoqvist et al visibility for Gibbsian spin-i systems 
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FIG. 88. Sjoqvist et al visibility for Gibbsian spin-1 systems 



FIG. 89. Sjoqvist et al visibility for Gibbsian spin-| systems 
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FIG. 90. Sjoqvist et al visibility for Gibbsian spin-2 systems 



FIG. 91. Sjoqvist et al visibility for Gibbsian spin-| systems 
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FIG. 92. Sjoqvist et al visibility for Gibbsian spin-3 systems 



FIG. 93. Sjoqvist et al visibility for Gibbsian spin-| systems 
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FIG. 94. Sjoqvist et al visibility for Gibbsian spin-4 systems 



FIG. 95. Sjoqvist et al visibility for Gibbsian spin-| systems 
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FIG. 96. Sjoqvist et al visibility for Gibbsian spin-5 systems 

VIII. COMPARISONS ACROSS GIBBSIAN A-LEVEL SYSTEMS OF SJOQVIST ET AL GEOMETRIC 

PHASES AND VISIBILITIES 



FIG. 97. Sjoqvist et al geometric phases for n-level Gibbsian systems (n = 2,..., 11) holding a = 1. At 9 = .25, the curve 
for n = 2 is dominant, followed in order by those for n = 4,6, 8,10, (all having positive values at 9 = .25) and n — 3, 5, 7, 9,11 
(all having negative values at 9 = .25), cf. Fig. |2l| 
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FIG. 98. Sjoqvist et al geometric phases for n-level Gibbsian systems (n = 2,... , 11) holding 9 fixed at In monotonically 
decreasing order are the curves for (the even) n = 2,4, 6, 8, 10, followed by those for (the odd) n = 3, 5, 7, 9, 11, cf. Fig. B2 




FIG. 100. Sjoqvist et al geometric phases for n-level Gibbsian systems (n = 2,..., 11) holding a = 2. In order of decreasing 
dominance at 9 = .25 are the curves for (the even) n = 2,4, 6, 8,10, followed by those for (the odd) n = 3, 5, 7,9,11, cf. Fig. M 
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IX. DIRECT COMPARISONS OF UHLMANN AND SJOQVIST ET AL RESULTS 


In Fig. 102, we plot for n = 6 the ratio of the Sjoqvist et al geometric phase (Fig. gl|) to the U hlma nn geom etric 
phase (Fig. ^), and in Fig. 11031 , the corresponding ratio for n = 11. In the next two figures (Figs. KI4 and 105 ), we 
show the ratios of the associated visibililities. 
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FIG. 103. Ratio 
Gibbsian systems 


of the Sjoqvist et al geometric phase (Fig. ^>6|) to the Uhlmann geometric phase (Fig. w for the eleven -level 


ratio 



1.5 


FIG. 104. Ratio of the Sjoqvist et al visibility (Fig. §) to the Uhlmann visibility (Fig. |l5|) for the sir-level Gibbsian systems 
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FIG. 105. 
systems 



X. SUMMARY 


We have reported an in-depth study here of two different methodologies for determining geometric phases for mixed 
states jl]0]. The analysis has been framed in terms of rotations (0(3)-orbits) of n-level Gibbsian systems, as originally 
proposed by Uhlmann ( 3 ). 

There are, of course, many features in these results. One of these is that the Sjoqvist et al geometric phases (sec. 1V|) 
appear to be less sensitive to the inverse temperature parameter (a) than do the Uhlmann geometric phases (sec. |nj). 

In particular, we have found (sec. IV) that as the number of levels of the Gibbsian density matrices increases (in the 
sample we have studied) from n = 2 to 11, the Uhlmann methodology often yields simple monotonic behavior. On 
the other hand, this seems to be largely absent using the Sjoqvist et al methodology (sec. VIII ), which appears among 
other things to be sensitive to the bosonic (odd n ) or fermionic (even n) character of the system under consideration, 
It is certainly important to note that if one sets the parameter a — given in ( 0 ) — equal to zero in the approach 
of Uhlmann (sec. [|), the resulting holonomy invariant (Tj) is simply equal to (—l) n+1 times the holonomy invariant of 
Sjoqvist et al [0, eq. (15)]. This would appear to help to explain the added complexity of the plots of the Uhlmann 
geometric phases (sec. |J) vis-a-vis those of the Sjoqvist et al plots (sec. ©• 

We have found compelling numerical evidence that for odd n the central (|~l|]-th) of the n eigenvalues (ordered 
in terms of absolute value) of the Uhlmann holonomy invariant (jTj) is always a real number (sec. |V C| ). (By the 
fundamental theorem of algebra, for odd n , one of the n eigenvalues must of course be real — since complex roots 
come in conjugate pairs — but clearly not necessarily the centrally located eigenvalue.) Also, the ([§], |~]|] )-entry 
of the Uhlmann holonomy invariant ( 0 ) is itself real, for odd n. In terms of absolute values, the diagonal entries of 
the invariant ( 0 ) appear to be always monotonically decreasing from the (upper left) (l,l)-entry to the (lower right) 
(n, n)-entry. (We note that the angular momentum operator J z — entering in our equations ( 0 ), ( 0 ) and (|l5l) — itself 
has declining [real] diagonal entries, that is, n — l/over2 ,..., — 


n— 1 ^ 
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